Electronic standing waves on the surface of the topological insulator Bi2Te3 
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A line defect on a metallic surface induces standing waves in the electronic local density of states 
(LDOS). Asymptotically far from the defect, the wave number of the LDOS oscillations at the 
Fermi energy is usually equal to the distance between nesting segments of the Fermi contour, and 
the envelope of the LDOS oscillations shows a power-law decay as moving away from the line 
defect. Here, we theoretically analyze the LDOS oscillations close to a line defect on the surface 
of the topological insulator Bi2Te3, and identify an important preasymptotic contribution with 
wave- number and decay characteristics markedly different from the asymptotic contributions. The 
calculated energy dependence of the wave number of the preasymptotic LDOS oscillations is in 
quantitative agreement with the result of a recent scanning tunneling microscopy experiment [Phys. 
Rev. Lett. 104, 016401 (2010)]. 



PACS numbers: 68.37.Ef, 73.20.-r, 73. 20. At 

I. INTRODUCTION 

Distinct surface-electronic properties, potentially rel- 
evant for spintronic applications, arise from the strong 
spin-orbit interaction in three-dimensional topological in- 
sulators (SDTIs^. Although the bulk electronic struc- 
ture of these materials resembles that of standard band 
insulators with electronic bands separated by an energy 
gap, the valence and conduction bands of the surface 
states form a conical dispersion and touch at the cen- 
ter of the surface Brillouin zone. These gapless surface 
states lack the standard twofold spin degeneracy, they 
are protected against backscattering, and the spin ori- 
entation of each plane-wave surface state is determined 
unambiguously by its momentum vector. 

In the past few years, surface-sensitive experimental 
techniques have been utilized to explore the remarkable 
properties of the surface electrons in 3DTIs. The lin- 
ear, Dirac-cone-like electronic dispersion and deviations 
from that were observed in various 3DTI materials using 
angle-resolved photoemission spectroscopjBHTl (ARPES), 
and the correlation between spin and momentum was 
demonstrated by the spin-resolved version of the same 
techniqu^. The role of electron scattering off point- 
like impurities and line defects on 3DTI surfaces, highly 
relevant for future attempts to design electronic devices 
based on these materials, has been studied via scanning 
tunneling microscopy (STM^Pm. In the vicinity of ob- 
stacles on the surface, characteristic standing wave pat- 
terns are formed due to the interference of initial and 
final scattering state^^^. These electronic standing waves 
contribute to the local density of states (local DOS, 
LDOS), therefore real-space mapping of them is possi- 
ble via STM. Theories describing the standing waves on 
3DTI surfaces have also been formulated recentljff^ltSil 

A line defect has translational symmetry in the di- 
rection it stretches along, hence the electronic standing 
waves in its vicinity are essentially one-dimensional (ID), 



(i.e., the LDOS varies only along the axis perpendicular 
to the line defect). This simple ID character of the in- 
duced LDOS pattern implies a relatively straightforward 
experimental and theoretical analysis of the effect, which 
serves as a strong motivation to consider such arrange- 
ments. A line defect arises naturally at the edge of a step 
formed by an extra crystal layer on the surface ,'^'^^'0^ 
hence this ID setup is accessible experimentally. 

Information on the electronic system can be extracted 
from the asymptotic decay exponent and wave number 
of LD OS o scillation around line defects. Theoretical 
j.ggy^^^i2ii8H20l indicate that the LDOS oscillations on the 
surface of a 3DTI, asymptotically far from a line defect 
and within the energy range of linear dispersion, decay 
with the distance x from the defect as x^'^/^. This de- 
cay exponent is in contrast with the ^ x^^^"^ decay seen 
in a standard two-dimensional electron gas,-^ and arises 
as a consequence of the absence of backscattering char- 
acteristic of surface electrons in 3DTIs. Recent STM 
data agrees with this prediction.!^ The wave vector of 
the asymptotic LDOS oscillations is usually equal to the 
distance between nesting segments of the constant-energy 
contour (CEC), which is the diameter of the Fermi circle 
in the above-mentioned case. This has been used in a 
recent experimeniP to confirm the linear dispersion and 
to infer the Fermi velocity on the Dirac cone in Bi2Te3. 

For energies well above the Dirac point, the topologi- 
cal surface conduction band of Bi2Te3 is subject to strong 
hexagonal warping. S TM d ata corresponding to this en- 
ergy range is availabl^^^^, however, the rather complex 
geometry of the dispersion relation has so far prevented 
an unambiguous theoretical interpretation of the obser- 
vations. In this work, we provide a theoretical investi- 
gation of LDOS oscillations created by a line defect on 
the surface of Bi2Te3. We describe the effect in an ex- 
act scattering-theory frameworlP^, yielding results that 
are not restricted to the spatial region asymptotically far 
from the defect, but hold also in the vicinity thereof 
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FIG. 1: (a,b) Surface-state conduction-band dispersion of 
Bi2Te3 along the FM (solid line) and TK (dashed line) di- 
rections of the surface Brillouin zone. (c,d) Surface-state 
conduction-band density of states of Bi2Te-i. For (a) and (c), 
Eq. ^ was used with parameter values^ = 2.55 eVA, 
A = 250 eVA^ and a,7 = 0. For (b) and (d), the param- 
eter values vo = 3.5 eVA, A = 150 eVA^, a = 21 A^ and 
7 — —19.5 eVA^ were used. Red crosses (blue points) rep- 
resent the measured dispersion along FM (TK) (data taken 
from Ref. Ilip . The zero of energy corresponds to the Dirac 
point of the spectrum. Energy intervals overlapping with the 
bulk valence band (BVB) and bulk conduction band (BCB) 
are also shown. 

This enables us to directly compare our results with ex- 
perimental data, the latter being usually taken close to 
the defect where features of the LDOS are most pro- 
nounced. In the energy range of strong hexagonal warp- 
ing, we identify a significant pre-asymptotic contribution 
to the LDOS oscillations, with wave number quantita- 
tively matching that of a recent experimentfSl. 

II. BAND-STRUCTURE PARAMETERS. 

In order to base our forthcoming calculations on an ac- 
curate surface-band dispersion, we first establish accurate 
values of the relevant band-structure parameters (defined 
below) of Bi2Te3. ARPES measurements'^ indicate that 
the surface bands of 3DTIs with the crystal structure of 
Bi2Te3 are subject to hexagonal warping, which can be 
described by the envelope-function HamiltoniarP^ 

H{k) = + VI, {k^cjy - kya^) + y - fc3 ) rr^ , (1) 

where = vo{l + ak^), and k± = i: iky. Here, 
(cr^, (7y, a^) is the vector of Pauli matrices represent- 
ing spin, vo,X,j and a are band-structure parameters, 
k — {kx'i 1 and kx and ky are momentum components 
along the FM and TK directions of the surface Bril- 
louin zone, respectively. For convenience, we performed a 
+TT /2 rotation around the z axis compared to the Hamil- 
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FIG. 2: (a) An incident wave from the a; > region (solid 
arrow) is partially reflected (dashed arrow) and transmitted 
(dotted arrow) at a line defect, e.g., an atomic terrace, on the 
surface of Bi2Te3. (b) Hexagonally warped constant-energy 
contour in reciprocal space, corresponding to energy E — 330 
meV above the Dirac point. Incident and reflected wave vec- 
tors from (a) are also shown. 

tonian in Ref. [HI Energy eigenvalues of H in Eq. ([T]) 

are 

e± (k) = 7k2 ± v'(t.kk)VA2fc2(fc2_3fc2)2 ^ (2) 

where + (— ) stands for conduction (valence) band. 
Higher-order term^2ll ^ can also be included in H . 

Satisfactory agreement between the spectrum in 
Eq. ^ and the ARPES spectra of Bi2Se3 surface stated 
can be obtained by neglecting the band-structure param- 
eters 7 and a. In Bi2Te3 however, the conduction-band 
dispersion measured along the FM direction, shown with 
red crosses in Figs.[l^,b, has a sub-linear segment, which 
can be theoretically reproduced by Eq. ([2| only if 7 and 
a are finite. We find that using the band-structure pa- 
rameter set Wo = 3.5 eVA, A = 150 eVA^, a = 21 A^ 
and 7 — —19.5 eVA^, the measured dispersion relations 
along the TK and FM directions and the surface density 
of states (Fig. lc,d of Ref. |TT1) are accurately described 
by H [Eq. ([l])] up to 335 meV above the Dirac point. 
We use these parameter values throughout this paper. 
Fig. (Tji and [TJd (Fig. ^ and [iji) compares the theoreti- 
cal surface-state dispersion (density of states) calculated 
with a parameter set used in an earlier worlP31, and the 
above parameter set that we found to be optimal, respec- 
tively. Further considerations used to find the optimal 
parameter set above are included in Appendix [A] 

III. MODEL. 

Our goal is to theoretically describe the oscillations in 
the surface-state conduction-band LDOS induced by a 
line defect, e.g., the edge of an atomic terracJsHll!^ on the 
surface of Bi2Te3. For the moment we assume that the 
defect forms a straight line that coincides with the y axis 
(see Fig. [2]). Following Ref. 18, we model the system 
with the Hamiltonian H + V, where effect of the defect 
is described via the potential V{x) = VbO(— x). 

Our analysis of the LDOS oscillations is based on exact 
energy eigenstates describing scattering of conduction- 
band electrons by the line defect (see Fig. [2|. Therefore 
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we first describe the scattering process of a plane wave 
energy eigenstate ^k,q{x,y)e^'^y = e^'^ye^'^^Xk,q incident 
from, say, the a: > side of the edge, with momentum 
components k = and q = ky, energy E, and spin 
wave function Xk,q- Scattering at the hne defect is elas- 
tic, hence the energy E is conserved. The momentum 
component q parallel to the defect is also conserved due 
to translational invariance in the y direction. The value 
of q determines the number of propagating waves at a 
given energy. For example, in Fig. [2]d, the number of 
propagating waves can be two (I and R) or four (Ii, I2, 
Ri, R2), depending on q. 

However, the incident plane wave can be scattered into 
coherent superpositions of three reflected and three trans- 
mitted partial waves, for the following reasons. On the 
X > side of the defect, the equation E = £+(fcr, q) has 
six complex solutions fc^.i, . . . fc^.e for given values of E 
and q, which follows from Eq. ([2| (the numerical method 
used to obtain these solutions is described in Appendix 
[b| . Three of the kr^p-s correspond to propagating waves 
moving away from the defect or evanescent modes. The 
associated wave functions ^k,.p,q = e*''^e*'^''-p'^x(fcr,pj <?) 
(p = 1,2,3) should be included in the Ansatz of the com- 
plete scattering state. The remaining three solutions fcr,p 
(p = 4, 5, 6) correspond to propagating waves towards the 
defect or diverging modes, hence they are disregarded. 

These arguments, together with their generalization to 
transmitted waves, imply that the a;-dependent compo- 
nent of the complete scattering wave function is 



where 
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1 otherwise. 



(4) 
(5) 



Here, the r-s and t-s are reflection and transmission coef- 
ficients, v±^k.q is the x-component of the group velocity of 
the plane wave with wave- vector components (fc, g), and 
the factors a and /3 ensure the unitarity of the scatter- 
ing matrix built up from the reflection and transmission 
coefficients. Evanescent modes are not subject to the 
unitarity requirement, hence we are allowed to make the 
above arbitrary choice a — (3 — 1 for partial waves with 
complex wave vectors. As the Hamiltonian is a third- 
order differential operator, partial waves at the two sides 
of the defect should be matched via boundary conditions 
ensuring their continuity as well as the continuity of their 
first and second derivatives. The scattering state ip'j,^^ of 
a plane wave incident from the left side (x < 0) of the 
defect can be described analogously. 



The LDOS at a given energy E and position x is ex- 
pressed with the exact scattering states as 



piE,x) 



where F^^ [F 
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is the wave- vector contour of waves that 
(i) are incident from the right [left] side of the line de- 
fect, and (ii) have energy E. Note that F^^ breaks up to 
disconnected pieces for strong hexagonal warping, e.g., 
in Fig. 2b, the points Ii and I2 belong to F^-* but Ri 

— ( R) 

and R2 do not. (The wave- vector contour F^ is shown 
in Fig. ^p, there it is formed as the union of the thick 
blue lines.) The infinitesimal line segment along r^'^-* 
is denoted by dn, and v(fc, q) is magnitude of the group 
velocity of the wave with momentum vector (fc, q). Using 

the exact scattering wave functions ip^g, we evaluate the 
integral in Eq. ^ numerically. To account for the in- 
evitable roughness of the line defect and to suppress noise 
due to the limited precision of the numerical integration, 
we average the LDOS oscillations Sp{x) over the angular 
range [—5°, 5°] of the line defect orientation with respect 
to the TK direction. 



IV. RESULTS. 

In Fig. [3^, we plot the numerically computed LDOS 
oscillations Sp{x) = p{x) — pq (black points) on the x > 
side of the defect, corresponding to energy E = 330 meV 
(3) and potential step height Vq = — 1 50 meV. Recent theo- 
ries using asymptotic analysi^i^'^oi g^ggest that the LDOS 
oscillations in the vicinity of a line defect on the sur- 
face of Bi2Te3 decay no faster than x~^^^. Motivated by 
this finding, we fit the function fi{x) = sin(2fcia; -|- 
Lpi)x^^^^ via fitting parameters Ai, ki, and tpi to the 
results (shown as red solid line). We also fit an exponen- 
tially decaying function f2{x) — A2 sm{2kfitx + ip2)e~^/^ 
via fitting parameters A2, fcfit, (P2 and L (blue dashed 
line). The two major features of our numerical result 
Sp{x) are as follows, (i) Comparison of the three curves 
suggests that the decay of dp{x) is better described by 
the exponentially decaying function f2ix) than by fi{x) 
having power-law decay (see Appendix [C] for further de- 
tails), (ii) The wave- number value obtained from fitting 
f2ix) is 2A:fit « 0.178 A'^. 

Expectations for the wave number of the LDOS oscilla- 
tions can be drawn from asymptotic analysiJ^^^. That 
suggests that the wave number of an electronic stand- 
ing wave at a given energy E is associated to wave vec- 
tors connecting nesting segments of the corresponding 
CEC"i4 2j - g^^ 2krM = 0.297 A"! or k^est = 0.126 
A~^ depicted in Fig. [3p. As the wave number 2A;fit char- 
acteristic of our data Sp{x) deviates significantly from 
2fcrM and fcnest, and its decay is exponential rather than 
power-law, we conclude that Sp{x) is dominated by a pre- 
asymptotic contribution in the considered spatial range. 
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FIG. 3: (a) Position-dependent contribution 5p{E, x) (solid line) to the LDOS at energy E — 330 meV in the vicinity of a line 
defect (black points). Red solid and blue dashed lines are fits of the functions /i and /2 (see text), respectively, to the data. 
po = 628 meV~^/im~^. (b) Constant-energy contour at the same energy, and relevant scattering wave vectors in reciprocal 
space. Thick blue (thin red) pieces of the contour correspond to left-moving (right-moving) plane waves. The union of the thick 
blue pieces form in Eq. (6]|. (c) Dominant wave number 2fcfit of 5p{E,x) (open circles), and characteristic wave numbers 
fcncst (red points), fcv (green diagonal crosses) and 2fcrM (blue crosses) of the constant-energy contour, as functions of energy 
E. (d) Magnitude of reflection amplitude \r\ = \rk.q'\, (e) magnitude of spinor overlap \x^x\ = Ixl q'X-k,q'\, (f) magnitude of 
the inverse of the parallel-to-defect group velocity component ii||(fc) (in units of 10~®s/m), and (g) the product of the above 
three quantitites (in units of 10~^s/m), as functions of perpendicular-to-defect wave number component k. 



In what follows, we argue that (i) the pre-asymptotic 
component of the LDOS is due to the interference of 
incoming and reflected partial waves, i.e., the role of 
evanescent and transmitted partial waves is negligible, 
and (ii) the appearance of the characteristic wave num- 
ber 2A:fit in the LDOS oscillations is due to the non- 
monotonic behavior of the parallel-to-defect group veloc- 
ity component along the CEC. To this end, we now con- 
sider only the interference contribution of incoming and 
reflected propagating waves to the LDOS [Eq. ^] on 
the right half plane p{E,x > 0), rewrite it as an integral 
over the perpendicular-to-defect wave-vector component 
fc, and drop the contributions from fc-regions where more 
than one propagating reflected partial wave is allowed 
(|fc| < fee), yielding 



Pr{E,x) 



dk 



{rk^q'xl.g'X-k^q'e ^^''"^ + c.c, 
\ME,k)\ 



(7) 

Here, q' = q'(k,E) is the unique positive solution of 
E = e+(fc, q) for a fixed k and E. The integrand without 
the exponential factor is related to the Fourier trans- 
form of Pr{x) . We plot the three factors determin- 
ing pr{x) — the magnitudes of the reflection coefficient 
|r^. g'l, the spinor overlap jx^, qiX-k,q'\i and the inverse of 
the parallel-to-defect group- velocity component v\\ (k, E), 
— as well as their product, in Fig. [Sji, e, f, and g, respec- 
tively. While Figs. [3ji and e show a featureless depen- 
dence on k, Fig. [sj^ reveals a peak in l/w|[(A:). The cor- 



responding local maximum point, which we denote with 
k^j , is very close to fcgt obtained from the numerical result 
in Fig. |3^. As a consequence of the peak in l/wj|(fc), a 
peak arises in the product of the three factors (Fig. [3^) 
as well. This analysis reveals that the characteristic wave 
number fcgt of the LDOS oscillation Sp{x) is determined, 
to a large extent, by the electronic dispersion relation via 
W|[ (k), and the details of the scattering process have little 
significance on its value. 

We have repeated the above analysis for various energy 
values in the range E S [145 meV, 475 meV] in order to 
compare the characteristic wave number 2fcfit of Sp{x) 
with experimental datgP^, and to confirm the correlation 
between the characteristic wave numbers obtained from 
the dispersion relation [kv{E)] and from the numerical 
LDOS calculation [kfit{E)]. We plot 2fcfit as the function 
of energy E in Fig. [Sj;. For low energy E < 170 meV, 
the hexagonal warping of the CEC is weak, and our re- 
sult shows fcflt ~ krM and a decay of Sp{x) oc a;~^/^ (not 



shown in the figures), in ag reement with the results of 
the asymptotic analysi^''^'-"'. Between 190 meV and 345 
meV above the Dirac point, our fc^t data in Fig. [3]: differs 
significantly from k-pAi , and the former shows good quan- 
titative agreement with the experimental values (shown 
in Fig. 4b of Ref. [TT]). Remarkably, 2ky{E), shown as 
green diagonal crosses in Fig. |3j:, is almost perfectly cor- 
related with 2fcfit, confirming the generality of the above 
proposition that the characteristic wave number of the 
LDOS oscillation is determined by the electronic disper- 
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No experimental data is available below 190 meV, 
whereas above 345 meV, the measured data (shown in 
Fig. 4b of Ref. [TT]) shows a pronounced kink that is 
not described by our model. A potential reason for that 
discrepancy is that the surface and bulk conduction elec- 
trons might be strongly hybridized in that high-energy 
range, making our surface-band model inappropriate to 
describe the corresponding standing-wave patterns. 

At high energy E > 200 meV, the nesting of CEC 
segments connected by fcnest in Fig. [3|d ma y also induce 
LDOS oscillations with wave number fcnost'^^'^- How- 
ever, in our model we find that such oscillations do not 
exist, due to the exact cancellation of contributions from 
reflected and transmitted waves incident from the x > 
and a; < regions, respectively (see Sec. |v|. 

The appearance of the LDOS contribution with wave 
number corresponding to the local maximum point of 
l/w|[(fc) is a generic feature, expected to be present in 
other electronic systems as well. We think that it plays 
a dominant role in Bi2Te3 because of the suppression of 
the other two Fourier components with wave numbers 
knest and 2fcrA/, due to the cancellation mechanism (see 
Sec. |v]) and the absence of backscattering, respectively. 



V. IRRELEVANCE OF TRANSMITTED WAVES 
TO THE LDOS OSCILLATION 

In principle, plane waves incident from the x < re- 
gion can contribute to the LDOS oscillations in the a; > 
region, provided that they are transmitted into at least 
two propagating channels on the x > side of the line 
defect. In this section, we show that such a contribu- 
tion is precisely balanced and canceled out in our model 
by the interfering reflected components of plane waves 
incident from the a; > side. In turn, this cancella- 
tion is responsible for the absence of LDOS oscillations 
with wave number A:,icst- Without the above cancella- 
tion mechanism, such oscillations would be expected to 
arise as fcncst connects nesting segments of the CEC (see 
Fig. 3b ). 

As we show now, our above statements follow from 
the unitary character of the scattering matrix describing 
the line defect. Following the notation used in Sec. HI 



consider electron plane waves with energy E and parallel- 
to-defect wave vector component q. We assume for con- 
creteness that for these given parameters E and q, there 
exists two incoming, and correspondingly, two outgoing 
plane waves on the a; > side, and one incoming and 
one outgoing wave on the a: < side. The corresponding 
perpendicular-to-defect wave- vector components are kn, 
ki2, koi, ko2 on the a; > side and ki and ko on the a; < 
side, respectively. In this example, the scattering matrix 
S{E, q) has the following structur^2lJ 



S 



11 



r 



12 



(8) 



A specific example is shown in Fig. where the points 
I\ O fcii and I1 O k\2 represent the incoming waves from 
the a; > region and i?i -H> fcoi and i?2 ^02 represent 
the reflected and transmitted waves in the a: > region. 
Note that the conclusions of this Section hold for different 
number of scattering channels as well. 

Each of the two electron waves incoming from the 
X > side is reflected into two propagating states with 
reflection amplitudes r'^,^ ip,p' — 1,2). The incident 
wave from x < is transmitted into two propagating 
states on the x > side with transmission amplitudes 
tp' {p' — 1,2). Straightforward calculation shows that 
the contribution of the transmitted waves to the LDOS 
oscillations in the x > region is accompanied by a con- 
tribution from the interference of the reflected waves, and 
the sum of these contributions is proportional to 



11' 21 



+ r T^r. 



12' 22 



11-2) 



t2 r: 



21 '"22 



(9) 

The first factor in Eq. ^ is the scalar product of the 
first and second rows of the scattering matrix describing 
the line defect, hence it vanishes because of the unitary 
character of the scattering matrix. 

This finding has the following remarkable consequence 
with respect to our calculated LDOS oscillations 5p{x). 
In the q intervals where two incoming and two outgoing 
waves exist (an example with a specific q is shown in 
Fig. (2)3, h o hi, I2 O fci2, Ri O fcoi and R2 O fco2) the 
wave number fco2 — fcoi approaches fcncst in a stationary 
fashion as q approaches its extremal value on the CEC, 
which, in principle, implies that the wave number fcnost 
is visible in the LDOS oscillations. In practice, however, 
the prefactor of the term oscillating with fco2 — fcoi, i-e., 
the first factor in Eq. ([9]), is zero for the whole q range 
with multiple outgoing waves. 

Note that the wave numbers fcip — fcop' (p,p' — 1,2), 
which appear in the LDOS oscillations due to interference 
between incoming and reflected waves, do approach fcnost 
as q approaches its extremal value on the CEC, but not 
in a stationary fashion. Consequently, these interference 
terms are also not able to promote fcnost to the dominant 
wave number of the LDOS oscillations. In summary, both 
the theoretical findings presented in this section and our 
numerical results shown above indicate that in the con- 
sidered parameter range, the characteristic wave number 
of the LDOS oscillation in the vicinity of the line defect 

is not fcnost- 



VI. DISCUSSION AND CONCLUSIONS 

A line defect on a metallic surface is usually mod- 
eled as a pote ntial stepP^EsHH or as a localized potential 
barriepISElIMl (e.g., Dirac-delta potential). In our work, 
we make the convenient, but arbitrary choice of modeling 
the defect as a potential step. It is important to note that 
the fine details of the LDOS results might in fact depend 
on the choice of the model of the defect (step vs. local- 
ized barrier) . However, our interpretation explaining our 
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FIG. 4: (a) Numerically obtained LDOS oscillations (points) 
at energy E = 330 meV for the case when the line defect 
is perpendicular to the VK direction, and the fit of f{x) = 
A8m{kx + (fi)/x^^^ with fitting parameters A, k and ip (solid 
red line), (b) Constant energy contour and the relevant wave 
numbers fcncst and 2krK- The wave number of the oscillation 
in (a) is given by fc^^t, as the oscillations with 2krK are fast 
decaying (oc x~^^^) due to the absence of backscattering. 



main result, i.e., the dependence of the standing wave's 
wave number kfit{E) on the Fermi energy E, is based on 
the momentum-dependence of the inverse velocity l/wy 
(see Figs. [3j;-g). The latter quantity is independent of 
the model describing the line defect, therefore our main 
conclusion is expected to hold even if (a) the height of the 
potential step at the line defect, denoted by Vb, is varied, 
or (b) a different model (e.g., Dirac-delta potential bar- 
rier) is used to represent the defect. To further support 
the expectation (a), we have numerically calculated the 
LDOS oscillations for various values of the potential step 
height Vq and found no qualitative change in the inferred 
kfit{E) data. 

Even though the line defect in the considered experi- 
ment Rcf. 11 was perpendicular to the TM direction of 
the surface Brillouin zonc^^, it is instructive, and regard- 
ing future experiments, potentially useful to consider the 
other high-symmetry case when the line defect is oriented 
perpendicular to the TK direction. In Fig. [4] we demon- 
strate that both the wave vector and the decay char- 
acteristics of the LDOS standing waves we obtain from 
our numerical technique are in complete correspon dence 
with the analytical results of asymptotic analysi^^^M!^ 
Namely, the wave number of the oscillation is given by the 
extremal (maximal) perpendicular-to-defect width k^^^ 
of the constant energy contour, whereas the decay goes 
as 5p{x) cx Xj'Jx. Oscillations with 2fcrK are not seen 
in Fig. presumably because they are fast decaying 
(oc x^^i^ due to the absence of backscattering^. 

In conclusion, we theoretically described pre- 
asymptotic electronic LDOS oscillations in the vicinity 
of a line defect on the surface of Bi2Te3, with wave num- 
ber and decay characteristics markedly different from the 
asymptotic ones. The calculated energy dependence of 
the characteristic wave number of the LDOS oscillations 
is in line with STM data. In a general context, our study 
highlights the importance of pre-asymptotic calculation 
of the surface-state LDOS oscillations in the analysis 



and interpretation of STM experiments. 

Note: While completing this manuscript, we became 
aware of a related workP^ on electronic standing waves on 
3DTI surfaces. Our results partially overlap with those 
in Ref. [211 Apart from various details of the model, 
the major distinctive features of our work are (i) the in- 
terpretation of the results in terms of the properties of 
the group velocity (Fig. [sji-g), and (ii) the quantitative 
agreement with the experimental results reported in Ref. 

EH 
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Appendix A: Band-structure parameters of the 
homogeneous system 

To base our calculation on an accurate surface-state 
dispersion relation, in Sec. |ll] we estimated the band- 
structure parameters of Bi2Te3 by comparing the the- 
oretical dispersion (see Eq. ([2])) and DOS with the ex- 
perimentally observed ARPES and STM spectra and the 
DOS derived from those. The four band-structure pa- 
rameters are vq, A, 7 and a. Here we outline the consid- 
erations we used for those estimates. 

The signs of 7 and a can be determined by considering 
the ARPES spectrum along the TM direction, shown as 
red crosses in Fig. [T|a, and b. Note that this cut of the 
dispersion relation corresponds to the function e^{kx, 0) 
in Eq. ([2]). The measured dispersion is linear for small 
wave number, its slope first becomes smaller as the wave 
number is increased, but then the slope increases again as 
wave number is increased further. This characteristic is 
naturally captured by a polynomial of the wave number 
with negative second-order coefficient and positive third- 
order coefficient. Since the third-order Taylor series of 
e+(fcj;,0) in kx around zero has the form 



.(fcj.,0) w vok^ + -ikl + av^kl 



(Al) 



we can conclude that 7 < and a > is required to 
describe the measured dispersion. The signs of the re- 
maining two parameters and A has no effect on the 
spectral properties, therefore we assign a positive sign to 
them. 

Having the signs of band-structure parameters estab- 
lished, we determined the values of the four parame- 
ters by systematic visual comparison of the experimen- 
tal dispersions (Figs. [I^,b), the DOS data obtained from 
ARPES and STM measurements (Figs. lc,d of Ref. ITT|) . 
and the corresponding theoretical curves. 
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FIG. 5: Comparison of the discrete Fourier transforms of (a) 
the computed LDOS osciUation Sp{x) shown in Fig. [3|i, (b) 
the function f2{x) = A2sm{2kfitX + (p2)e~^'''" , where A2, fcflt, 
ip2, L are obtained from fitting /2 to 5p, (c) yl2 sin(2fcflta:: + 
ip2)x~'^^^, and (d) A2 sin(2kfitX + <p2)x' 



-1/2 



Fig. [3^, the exponentially decaying /2 provides a better 
fit, hinting that the decay characteristics is clos er to an 
exponential than to a power-law predicted earliei^i^ElHSol 
However, our fitting procedure is not conclusive, as the 
exponentially decaying function /2 has an extra fitting 
parameter, the length scale L of the decay. 

In order to investigate the decay characteristics fur- 
ther, here we provide the discrete Fourier transform (FT) 
(Fig. [5^) of the real-space data in Fig. [3^, and compare 
that to discrete Fourier transformed oscillations that de- 
cay in an exponential (Fig. [5|d) or power- law (Fig[5j:,d) 
fashion. The discrete Fourier transformation is carried 
out after symmetrization of the real-space data, i.e., after 
mapping the real-space data set f{xi) {i = 0, . . . , N — 1) 
to /sym.j (j = 0, . . . , 2N — 1) via the definition 



fixj) if < j < AT - 1, 
f{x2N-i-j) iiN<j<2N-l. 



(CI) 



Appendix B: Plane-wave states 

We label the electronic plane-wave states by their en- 
ergy E and the parallel-to-defect wave-number compo- 
nent q, which are conserved in the scattering process on 
the line defect (for details see the main text). Here we 
review a numerical method to obtain these states, which 
is necessary to solve the scattering problem at the edge 
step. For a given E and q there are six solutions of lon- 
gitudinal wave vector kj-^p (p = 1 . . . 6) which satisfies the 

characteristic equation det H{kr,q) — IE =0, where I 
is the 2x2 identity matrix and H is the Hamiltonian 
defined by Eq. (1) in the main text. Complex roots fc^.p 
of the characteristic polynomial 



det 



H-Ie] = Y,a{E,q)k'; 



(Bl) 



fe=0 



are equal to eigenvalues of the companion matrix:-^ of 
this polynomial, hence we find the roots k^.p by numer- 
ically diagonalizing the companion matrix. Then, the 
spinor component x of the corresponding plane wave 
<i>fc^ p,q(x, y) = e**^e*'^'' p'^x(^r,p, 9) can be numerically 
computed from the eigenvalues problem 



for p — 1^ 



H{kr^p)x{kr,p,q) = Ex{kr.p,q) 
,6. 



(B2) 



Appendix C: Fourier analysis of Sp{x) 



In Sec. |IV[ we present numerical results for the defect- 
induced spatial modulation of the LDOS, 6p{x). In order 
to develop an understanding of the decay characteristic 
of the LDOS oscillations, we fit a power-law decaying 
{6p cx x~^/'^) as well as an exponentially decaying func- 
tion, /i and /2, respectively, to our data. According to 



This symmetrization ensures that the Fourier transform 
will be real valued in the large N limit. 

The subplots of Fig. [s] show the FT of (a) our nu- 
merical results Sp{x) shown in Fig. [3^, (b) the function 
f2{x) = A2 sin(2fcfita; -|- ip2)e''''/^, where A2, fcfit, V2, L 
are obtained from fitting /2 to 5p, (c) A2 sin(2/cfita; -I- 
ip2)x~^^'^, and (d) A2 sin(2fcfita; -I- ip2)x^^^'^. The data in 
Fig. [5];,d shows, in accordance with the known analyti- 
cal formula describing the Fourier transform of power-law 
decaying sinusoidal oscillationJiS!^ that the FT develops 
non-analytical behavior at the characteristic wave num- 
ber 2kfif In contrast, our data set in Fig. [5^ shows 
no such non-analytical behavior, similarly to the FT of 
the exponentially decaying oscillation in Fig. [Sfj. This 
observation, although still not conclusive, further sup- 
ports the possibility that the LDOS oscillations in the 
vicinity of the line defect do not follow a power-law de- 
cay, and perhaps are closer to an exponential. Further 
analyti cal stud ies, e.g., the extension of the asymptotic 
analysiji^llllini to the pre-asymptotic spatial region might 
help settling this open issue. 

In Sec. |IV[ we argued that the pre-asymptotic com- 
ponent of the LDOS is due to the interference of in- 
coming and reflected partial waves i.e., that the role of 
evanescent and transmitted partial waves is negligible. 
To strengthen that point further, we plot the quantity 



Pr{2k) 



Re \rk,q'Xk,q'X-k.,q' 

2iih\v\\{E,k)\ 



(C2) 
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[le 



as a function of k [for definitions, see around Eq. 
as a dashed blue line in Fig. [sja,. Note that Pr{2k) is t 
Fourier transform of the symmetrized Pr{x) of Eq. ([7|). 
Figure[5^ further supports the interpretation that the dip 
around the characteristic wave number 2fcfit of the LDOS 
oscillation forms as a result of interference of incoming 
and reflected partial waves. 
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